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A non zero, surprisingly large value of the third mixing angle 613 has been 
measured in reactor neutrino experiments. Currently the most precise measure- 
ment of sin^ 2013 has been performed by the Daya Bay experiment sin^ 2S13 = 
0.089±0.010(stat.)±0.005{syst.) (7.7 a significance of sin^ 26*13 > ),the RENO ex- 
periment has measured the value sin'^ 2(9i3 = 0.113±0.013(stat.)±0.019{syst.) (4.9 
a significance) and the Double Chooz experiment sin'^ 2613 = 0.109±0.030(stat.)± 
0.025(syst.) (2.9 cr significance). These results are extremely important for future 
searches of violation of combined CP parity in lepton sector of the Standard model. 

1 Neutrino mixing and oscillations 

Neutrino flavor eigenstates Vf originating in weak decays together with charged 
leptons of three known flavors f = e, ji^T are super-positions of three mass eigen- 
states i^i i — 1,2,3 with masses toi, TO2, ma. The mixing is described by 3a;3 unitary 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix U with elements 

Using this mixing matrix one can write flavor eigenstates as super-positions of mass 
eigenstates as follows: 



The CPT invariance implies the equivalence of masses of neutrinos and anti- 
neutrinos and relates the mixing matrices for neutrinos (U) and anti- neutrinos (U): 

The mixing of anti-neutrinos flavor and mass eigenstates therefore obeys the rela- 
tion: 



Experimentally firmly established phenomenon of oscillations of neutrino flavors 
can be explained by non-diagonal mixing matrix and different masses of neutrinos 
mass eigenstates. In general the 3x3 unitary matrix contains 9 free parameters. For 
Dirac neutrinos one can reduce the number of physical parameters to 4 by re-phasing 
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five out of six iepton fields. Canonical representation of the mixing matrix is the 
ordered product of three rotations with angles 612, 9i3, 623 and one CP violating 
phase 5. In case of Majorana origin of neutrinos the mixing matrix contains another 
two so called Majorana phases ai,2- Oscillations of neutrino flavors do not depend 
on Majorana phases; the relevant for oscillations part of the mixing matrix is: 

/ 1 \ / ci3 siae*'^ \ / C12 S12 \ 

C//* = C23 S23 1 -S12 C12 

V -S23 C23 / V -siae-*^ C13 / V « 1/ 
where Sij and Cy denote sin 6ij and cos 6ij respectively. 

With an approximation of Am|]^ = to| — m\ = m\ — m\ = Am|2 the survival 
probability for electron anti-neutrino with energy E at the distance x from the 
source is given by the formula: 

P... (I) - 1 - ^i"^ 'i-' (^) - (fe) sin^ (2.„) (^) 

Values of two mass square differences and mixing angles 612 and ^23 have been 
measured in following experiments: 

• experiments with atmospheric and accelerator neutrinos: 

|Am§J = |m§ - m^l = |Am§2l = |m§ ~ml\^ (48.2 meVf and 6*23 = 45°. 

• experiments with reactor and Sun neutrinos: 
Am^i =ml-ml'=^ (8.7 meVf and 6*12 = 34°. 

Because |Am§]^| = 30 • Arriji there are two oscillations with 30 times different 
distances to the 1st oscillation minimum List^^„: 

|Am|i| = |Am§2| ii'*„„„ = 0.5km/MeV = 500km,/GeV 
Aml^ ii=.t™n = l5km/MeV = 15000km/GeV 

The aim of three reactor neutrino experiments Daya Bay [2], Double Chooz 
[3] and RENO [3] was to measure not observed before disappearance of reactor 
electron anti-neutrinos at distances of 0.5 km/MeV. Such observation would be a 
discovery of non zero value of the mixing angle ^13. An alternative way to measure 
the value of ^13 is via the observation of the appearance of electron (anti-) neutrinos 
in accelerator muon (anti-)neutrinos at distances 500 km/GeV. 

2 Detection of reactor anti-neutrinos 

Nuclear power reactors are very powerful sources of electron anti-neutrinos. Fission 
products are neutron rich isotopes that are transformed to stable nuclei via series 
of [3~ decays. In average « 6 electron anti-neutrinos follow each fission producing 
the flux of « 2 • 10^° of Pe per second per 1 GWth of reactor thermal power. Energy 
spectrum of anti-neutrinos decreases with the energy and extends up to w 10 MeV. 
Anti- neutrinos with energies E^, > 1.8 MeV can be detected by charged current 
interactions (so called inverse beta decay IBD): 

+ P ^ + n 
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on target protons of liquid scintillator. Detected anti-neutrinos energy spectrum 
is the product of original anti-neutrino flux and neutrino energy dependent IBD 
cross section. Detected spectrum has a typical shape with a maximum at E^, w 
4 MeV. 

2.1 Detection of positrons 

Positron losses its kinetic energy by ionization and annihilates with an electron 
producing pair of me — 511 keV gammas. This so called prompt signal Ep 

Ep =Te+ 2me 

is related to neutrino energy as follows: 

E,, = Ep + 0.8 MeV 

2.2 Detection of neutrons 

Kinetic energy of IBD neutrons is ~ 4 keV and 40 keV for neutrinos with energies 
4 MeV and 10 MeV respectively. Neutrons are thermalized by elastic collisions 
with Hydrogen and Carbon nuclei of liquid scintillator and then captured. In 
order to distinguish neutron signal from backgrounds, liquid scintillator is doped by 
Gadolinium. Two isotopes 155 and 157 of natural Gd have huge thermal neutrons 
capture cross sections of 61 and 254 k barns respectively. Neutron capture on these 
Gd isotopes is followed by the emission of several gammas with total energies of 
8.54 MeV and 7.94 MeV respectively. Because the natural abundances of ^^^Gd 
and ^^"^Gd are high (14.8% and 15.7%) the addition of a small amount of 0.1% of 
Gd is enough to ensure a capture of more than 80% of neutrons on Gd (and less 
then 20% on H of liquid scintillator). Time difference At between prompt signal 
from positron and delayed signal from neutron follows an exponential distribution 
g-'^*/ i'^) with the value of capture time (t) = 29/is inversely proportional to Gd 
concentration. 

3 Reactor neutrino experiments 

Number of detected events is proportional to the product of the detector mass and 
power of nuclear reactors. In addition optimal sensitivity requires to place far de- 
tector(s) at the distance close to 2 km where the maximal disappearance is expected 
for neutrinos around peak energy of 4 MeV. Correlated systematic errors are can- 
celled and uncorrelated reactor uncertainties are minimized in a configuration with 
identical near detector (s) [7] placed closer to reactors where oscillations effects are 
small. The background can be minimized by large overburden. The near-far ar- 
rangement of Daya Bay anti- neutrino detectors (ADs) , is shown on Fig. [1] relevant 
parameters for all the three experiments are summarized in the following table: 





Reactor Power 


Detector 


mass [t] 


Baseline 


[m] 


Overburden 


[mwe] 


Experiment 




Near 


Far 


Near 


Far 


Near 


Far 


Daya Bay 


17.4 


40, 20* 


60* 


470, 576 


1650 


250, 265 860 


Double Chooz 


8.5 




8.2 


400 


1050 


120 


300 


RENO 


16.5 


16.5 


16.5 


409 


1444 


120 


450 
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Table shows the configurations at the time of the PIC2012 conference. 
* Since October 2012 the Daya Bay detector has been completed to 40+40 t of near 
detectors and 80 t of far detectors. 

** The Double Chooz experiment currently operates with far detector only. 

L3 



• L4 
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Daya Bay NPP 



Figure 1. Layout of the Daya Bay experiment. The dots represent reactor cores, labelled as Dl, 
D2, LI, L2, L3 and L4. Six anti-neutrino detectors (ADs) were installed in three experimental 
halls {EHs).The figure is taken from [2]. 



3.1 Anti-neutrino detectors 

Each of three near and three far Daya Bay anti-neutrino detectors jS] is the stainless 
steel vessel with three nested cylindrical volumes separated by two concentric acrylic 
vessels (see Figure [5]). The innermost volume serves as the target and is filled with 
20 t of Gd loaded hquid scintillator (Gd-LS). In order to fully contain gammas 
from positron annihilation and Gd de-excitation the central volume is surrounded 
by gamma catcher of 21t of Gd free liquid scintillator (LS) in outer acrylic vessel. 
Outermost volume of 37 t of mineral oil serves as shielding of inner volumes from 
radiation originating from the walls of stainless steel vessel and photomultipliers. To 
increase the light yield and improve the homogeneity refiector panels are installed at 
the top and bottom of the outer acrylic vessel. Scintillating light is detected by 192 
8-inch PMTs placed along the circumference of the stainless steel vessel. The Daya 
Bay detectors are submerged in instrumented water pools such that at least 2.5 m of 
water is surrounding each detector. Purified water passively shields detectors from 
the ambient and cosmic rays induced radioactivity and actively detects Cherenkov 
light produced by cosmic muons. In addition, an array of resistive plate chamber 
(RPC) modules is placed on top of the pools for increased background rejection. 
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Target mass of the Double Chooz detector (see Figure [5]) is 8.5 t (10.3 m'^) of 
Gd-LS, surrounded by 22 m'^ of LS gamma catcher and by 110 m^ of mineral oil 
buffer. The detector is instrumented by 390 10-inch PMTs. The Double Chooz 
detector is submerged in 90 m^ of liquid scintillator active inner veto and it is 
covered and surrounded by 15 cm thick shield of demagnetized steel. From top it 
is covered by two active outer veto detectors made out of scintillator strips. 

The RENO experiment operates one near and one far detector each with the 
target mass of 16 t surrounded by 30 t of gamma catcher and 55 t of mineral oil. 
The detectors are instrumented with 354 10-inch PMTs and submerged in 350 t of 
active purified water veto system. 




Figure 2. Schematic diagram of the Daya Bay (left picture) and Double Chooz (right picture) 
detectors. Pictures are taken from Ref. |2I5) . 



3.2 Detector calibration 

Detectors of all three experiments are equipped with calibration systems. Daya Bay 
detectors have three automated calibration units ACU A, B and C (see Figure [U 
that insert calibration sources along vertical direction in the centre of the detector, 
close to the edge of target volume and in the middle of gamma catcher; Double 
Chooz calibration sources are lowered from a glove box (see Figure[2]) at the detector 
top; RENO detectors are using three dimensional and one dimensional calibration 
system. 

Double Chooz and RENO detectors are calibrated by gammas from ^'^^Cs ra- 
dioactive source: 

"^Cs ^ i3^™Ba + e- + i>e 

"^'"Ba ^13^ Ba + 7(661.7 keV) 
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The calibration of all three experiments near the IBD threshold is done with 
^^Ge that decays by electron capture to /3+ emitter ^^Ga following the decay 
scheme: 

e" Ge -^^^Ga + 

^^Ga Zn + e+ + 

detectors response to pair of 511 keV gammas from positron annihilation in the 
source capsule is reconstructed. 

Near the peak of detected anti-neutrino energy spectrum detectors of all three 
experiments are calibrated with ^°Co radioactive source: 

^^^Co -> ^^Ni* + e- + Pe 

60^-* ^60 ^ ■ _^ ^(1173,2 keV) + 7(1332.5 keV) 

Neutrons from ^^^Am—^^C source (Daya Bay) and from spontaneous fissions of 
^^^Cf (Double Chooz and RENO) are used to calibrate detectors response to gam- 
mas from neutron capture on H and Gd. The detector's response is also monitored 
by neutrons produced by cosmic rays. 

In addition the PMT timing, gain and relative quantum efficiencies are cali- 
brated using Light-Emitting Diodes (LED) light (Daya Bay and Double Chooz) 
and Laser light (RENO). 

Energy resolution of Daya Bay detectors has been measured to satisfy the re- 
lation a/E = 7.5%/y/E[MeV] + 0.9%; energy resolution of RENO follows the 
relation a/E = 5.9%/ y^E[MeV] + 1.1%. Values of statistical terms of 7.5% and 
5.9% are in good agreement with the yield of 160 photoelectrons/MeV and 250 
photoelectrons/MeV measured by Daya Bay and RENO detectors respectively. 



3.3 IBD candidates selection 

Daya Bay. Minimum value of the IBD signal at the threshold is 1.022 MeV. Due 
to energy resolution and possible inhomogeneity of the response, prompt signal 
is selected in the region 0.7 MeV < Ep < 12 MeV. Delayed signal 6 MeV < 
Ed < 12 MeV from the neutron capture on Gd is searched in time window At = 
td — tp = 1 -r- 200 fxs following the prompt signal. In addition there should be 
no signal exceeding 0.7 MeV in 200/is time windows preceding and following the 
prompt signal and the delayed signal respectively (so called multiplicity cut). 

In order to suppress cosmogenic background Daya Bay detectors are vetoed for 
0.6 ms, 1 ms and 1 s after so called pool muon, non-showering AD muon and 
showering AD muon signals respectively. 



Selection criteria for all three experiments are compared in following table. : 



Detector 


Ep 


Ed 


td tp 


Multiplicity 


Veto 




MeV 


MeV 




jJLS 


ms 


Daya Bay 


(0.7; 12) 


(6; 12) 


(1; 200) 


200 / 200 


0.6 / 1 / 1000 


Double Chooz 


(0.7; 12.2) 


(6; 12) 


(2; 100) 


100 / 400 


500* 


RENO 


(0.7; 12) 


(6; 12) 


(2; 100) 


100 / - 


1 / 10** 



Two values for multiplicity cut denotes the time window preceding the prompt 
signal and following the delayed signal respectively. 
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* Double Chooz detector is vetoed for 500 ms after cosmic muon signal >600 MeV. 
** RENO detectors are vetoed for 1 ms after AD muon signal larger than 70 M eV 
(or >20 MeV if that coincides with large water pool signal) and for 10 ms after 
cosmic muon signal >1.5 GeV. 

Total number of IBD candidates collected by each experiment is summarized in 
the following table: 



Detector 


Near 


Far 


Day a Bay 
Double Chooz 
RENO 


205 308 
154 088 


28 909 
8 249 
17 102 



The correlation of daily IBD rates per detector and reactor powers for Daya 
Bay and RENO are shown on Figure [3l 
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Figure 3. The correlation of daily IBD rates and reactor powers. Left plot shows Daya Bay daily 
IBD rates per detector in three experimental halls (near EHl and EH2 and far EH3). Right plot 
shows RENO daily IBD rates measured in near and far detector. Figures are taken from Ref. I2I3I . 



3.4 Systematic uncertainties 

The values of systematic uncertainties for all three experiments are summarized in 
the Table [B 



Table 1. Comparison of systematic uncertainties. All values are in %. 



Detector 


Detec 
Corr 


;or related 
Uncorr 


React 
Corr 


or related 
Uncorr 


Back 
Si 

Near 


ground 
'inal 

Far 


Daya Bay 
RENO 

Double Chooz 


1.9 
1.5 


0.2 
0.2 

1.0 


3.0 
2.0 


0.8 
0.9 

1.7 


1.9 ± 0.2 
2.8 ± 0.8 


4.7 ± 0.35 

5.8 ± 1.1 
5.5 ± 1.6 



For Daya Bay and RENO experiments with near-far configuration of detectors 
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only uncorrelated systematic errors influence the result. Dominant contributions 
to detector related uncorrelated uncertainties are uncertainties of delayed energy 
cut efficiency and of fraction of neutron captures on Gd. Due to near-far detector 
configurations, influences of uncorrelated reactor based uncertainties to measured 
ratio of far to near IBD rates are highly reduced (by a factor of « 20 for the Daya 
Bay). 

3. 5 Backgrounds 

The signal of IBD events can be mimic by three major types of background events. 

• Accidental background are coincidences of prompt and neutron-like signals 
not caused by IBD events. This background is estimated using data and it is 
concentrated at low energies. 

• Fast neutrons background. Energetic neutrons created in photonuclear inter- 
actions of cosmic muons could mimic IBD by recoiling off a proton before being 
captured on Gd. The number of background events can be estimated by the 
extrapolation of the background shape at energies above 12 MeV. 

• Very neutron rich exotic nuclei created in interactions of cosmic muons. These 
are mostly ^Li and less frequently also |-ffe. Large fraction of ^Li and |7?e 
/3~ decays is accompanied by free neutron: 

iLi -¥ iBe + n + e- + i7e Q ^ 13.607 MeV 
iBe 2a 

iHe^lLi + n + e- +i7e Q ^ 10.651 MeV 

Due to large Q value f3^ — n decays of these nuclei mimic IBD candidates 
spanning the whole spectrum and because of their long half-lives of 178 and 
119 ms detectors shall be vetoed for long times after showering cosmic muon 
signals. 

The values of background to signal ratios and their uncertainties are shown 
in the last column of the Table [TJ In all three experiments the uncertainty of 
background to signal ratio is dominated by Li/He background; it is lowest for the 
Daya Bay due to overburden and strict veto cut (1 s) after muon interacting in the 
volume of anti-neutrino detectors. 

4 Results 

Daya Bay 

The De rate in the far hall was predicted with a weighted combination of the two near 
hall measurements assuming no oscillation. A ratio of the measured to expected 
rate is defined as 

^^Mf _ Mf 

Nf aMa + l3Mb' 
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Figure 4. Ratio of measured versus expected signals in six Daya Bay detectors (left plot 2 ) and 
two RENO detectors (right plot [3]). The oscillation survival probability at the best-fit value is 
given by the smooth curve. The value versus sin'^ 2Si3 are shown in the inset. Figures are 
taken from Ref. I2I3| . 



where N f and Mf are the predicted and measured rates in the far hall (sum of AD 
4-6), Ma and Mi, are the measured, background-subtracted IBD rates in detectors 
placed in EHl and EH2, respectively. The values for weights a and (3 were domi- 
nated by the baselines, and only slightly dependant on the integrated flux of each 
core. The ratio observed at the far hall was: 



R = 0.944 ± 0.007(stat.) ± 0.003(syst.), 



where the statistical (systematic) uncertainties were obtained by propagating sta- 
tistical (uncorrelated systematic) uncertainties in the measured IBD counts in the 
three halls. 

The value of sin^ 20i3 was determined by minimizing a functional constructed 
with pull terms accounting for the correlation of the systematic errors, see details 
in [2]. The best-fit value is 



sin^ 26*13 0.089 ± 0.010(stat.) ± 0.005(syst.), 



with a x^/NDF of 3.4/4. All best estimates of pull parameters are within its one 
standard deviation based on the corresponding systematic uncertainties. The no- 
oscillation hypothesis is excluded at 7.7 standard deviations. 
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4.2 RENO 

Rate based analysis of data taken by RENO experiment yielded following results 
for ratio of measured to expected rate in the far detector: 

Rfar = 0.920 ± 0.009(stat.) ± 0.014(syst.), 

and for the value of sin^ 2^13: 

sin^ 26*13 = 0.113 ± 0.013(stat.) ± 0.019(syst.), 

with the significance of 4.9 standard deviations. 

The results of Daya Bay and RENO experiments are illustrated on Figure ID 
Distortion of the observed ratio of far and near D,, spectra (see Figure [5]) provides 
further evidence of neutrino oscillations. 




Prompt Energy (MeV) ' Prompt energy [MeV] 



Figure 5. Left: Daya Bay. Top: Measured prompt energy spectrum of the far detectors compared 
with the no-oscillation prediction based on the measurements of the near detectors. Bottom: The 
ratio of measured and predicted no-oscillation spectra. The solid curve is the expected ratio with 
oscillations, calculated as a function of neutrino energy assuming sin^ 2^13 — 0.089 obtained from the 
rate-based analysis. The dashed line is the no-oscillation prediction. Right: R-ENO. Top: Measured 
prompt energy spectrum of the far detector compared with the no-oscillation prediction based on the 
measurements of the near detector. Bottom: The ratio of measured and predicted no-oscillation spectra. 
The dashed line is the no-oscillation prediction. Figures are takeir from Rcf. I2I3I . 



4-.3 Double Chooz 

In Double Chooz experiment running without near detector, the data from far 
detector are compared to calculated (and normalized to Bugey4 data) no-oscillation 
anti- neutrino spectra. The rate only analysis result is: 
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sin^ 26*13 = 0.170 ± 0.035(stat.) ± 0.040(syst.). 

Double Chooz collaboration reported also combined rate and shape analysis [5] 
result illustrated on Figure |6l 

sin^ 26*13 = 0-109 ± 0.030(stat.) ± 0.025(syst.). 




Energy [MeV] 

Figure 6. Double Chooz. Top: Measured prompt energy spectrum superimposed on the expected 
prompt energy spectrum, including backgrounds (green region), for the no-oscillation (blue dotted curve) 
and best-fit (red solid curve) at sin^ 2^13 — 0.109. Inset: stacked spectra of backgrounds. Middle: Ratio 
of data and no-oscillation prediction. Bottom: differences between data and no-oscillation prediction 
(data points), and differences between best fit prediction and no-oscillation prediction (red curve). The 
orange band represents the systematic uncertainties on the best-fit prediction. Figure is adopted from 
Ref. [5]. 



5 Importance of the result 

Non zero value of ^13 would allow for violation of CP and T symmetry in lepton 
sector of the Standard Model of elementary particles. In particular following CP 
and T violating differences of appearance probabilities can be measured {f ^ g)'- 




12 



R. Leitner 



Due to CPT invariance, they shall have the same values: 



For three neutrino flavours the variable A fg has the same absolute value for all 
f g that is given by formula: 



Using current measurement of 613 and other neutrino mixing parameters one can 
estimate the effect of CP violation to the difference of P^^^^^^ — Pu^^u, measured 
at the first appearance maximum [L/ E = 500 km/GeV) to be « 3% • sin (5 . 

6 Summary 

A non zero, surprisingly large value of the third mixing angle ^13 has been measured 
in 2012. The result is extremely important as it is opening future searches for 
violation of CP in lepton sector. 

After 2011 hints for non zero value of 6'i3 from accelerator experiments |8|l9j and 
Double Chooz [3] it is important that today we have consistent results from three 
different experiments: 

• precise measurement of ^13 by Daya Bay [2] experiment with the significance of 
7.7 sigma reported at Neutrino 2012; the discovery of non zero value with sig- 
nificance exceeding 5 sigma was announced in March and published in Ref . [1] . 

• 5 sigma observation of non-zero value of ^13 by RENO announced in April and 
published in Ref. [3] 

• 3 sigma indication of non-zero value of ^13 by Double Chooz [5] experiment 
reported at Neutrino2012. 

With more data in near future one can expect improved results with reduced 
statistical errors and systematic uncertainties. More detailed analysis will include 
also the shape of detected energy spectra. 

The three experiments collected several hundred thousands of anti-neutrino in- 
teractions at different distances from the reactor cores and new interesting analyses 
will be performed using such unique set of data. 
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